Trapped Two-Dimensional Fermi Gases with Population Imbalance 
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We study population imbalanced Fermi mixtures under quasi-two-dimensional confinement at zero 
temperature. Using mean-field theory and the local density approximation, we study the ground 
state configuration throughout the BEC-BCS crossover. We find the trapped system to be either 
fully normal (throughout the trap) or to consist of a superfluid core phase surrounded by a normal 
shell, which is itself either fully or partially polarized. Upon changing the trap imbalance, the 
entire trap configuration may undergo continuous or discontinuous transitions between the different 
ground states. We establish the phase diagrams for trapped situations both with and without a 
surface tension. The latter is found to strongly favor the fully polarized normal phase. We also 
study the collective excitations of the system and find discontinuities of the frequency spectra at 
the discontinuous transition point. There also appear modes with oscillation frequencies below 
the trapping frequency. Finally, we argue that thermal equilibration throughout the trap will be 
considerably slowed down at low temperatures when a superfluid phase is present. 
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I. INTRODUCTION 

The recent experimental progress on manipulating and 
cooling fermions has allowed the experimental realization 
and study of imbalanced fermion gases [E S S IE IE 13 • 
An important experimental tool is the trapping of ultra- 
cold atoms in optical lattices using which various solid- 
state systems may be simulated, while the dimensionality 
of the system may be reduced to 1 or 2 dimensions [8| . 
However, perhaps the most extraordinary experimental 
possibility is the tunability of the interparticle interac- 
tions using Feshbach resonances S E3- By adjusting 
the scattering length one is able to study the full BEC- 
BCS crossover in which weak interparticle attractions 
give rise to a Bardeen-Cooper-Schrieffer (BCS) type cou- 
pling while strong interactions lead to the formation of 
bound pairs which then undergo Bose-Einstein conden- 
sation (BEC). 

Understanding the ground state properties of the three- 
dimensional imbalanced three dimensional Fermi system 
is the subject of several recent theoretical works [ll|, EH 
El 0, El El El El EISqI- Zero-temperature Monte 
Carlo (MC) simulations renormalization group (RG) 
calculations [12l ] and a model which incorporates beyond- 
mean-ficld and nonzero temperature fluctuations [20| 
were all found to be in reasonable agreement with the 
MIT experiments H 0, H H 0] ■ The Rice group's exper- 
iments P, 0] , on the other hand, appear to be incompat- 
ible with the local-density approximation (LDA). Tak- 
ing into account surface effects in addition to the LDA 
is believed to be the key to understanding this disagree- 
ment d El) EH EH > alt houg h nonequilibrium phenomena 
may also play a role [2l|, |22| . Also recently, both theoreti- 
cal and experimental studies of collective excitations with 
imbalanced fermion gases were performed [23l. l24l l25l |2(| . 



The physics of two-dimensional ultracold gases turns 
out to be significantly different from that of three di- 
mensions. In practice, a quasi-2D system can be real- 
ized using optical lattices. This has already been done 
with balanced Fermi mixtures [1, [27| . To date, the- 
oretical studies of 2D imbalanced fermion gases have 
mainly focussed on homogeneous systems and were 
concerned with the effects of finite temperature 12811 . 
the Berezinskii-Kosterlitz-Thouless (BKT) transition [23, 
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we explore trapped, imbalanced 2D 
fermion gases at zero temperature usin g th e mean-field 
theory for homogeneous phases of Ref. [3y. We present 
phase diagrams which depend only on an interaction pa- 
rameter and the population imbalance. Three relevant 
trap configurations are identified: a fully normal config- 
uration (FN) and two configurations consisting of a su- 
perfluid (SF) core surrounded by a fully polarized normal 
(FPN) or partially polarized normal (PPN) shell, respec- 
tively. The transition between these three configurations 
can be either continuous or discontinuous. We also study 
the influence of a nonzero N-SF surface tension on the 
phase diagrams and find that it strongly suppresses the 
PPN phase. Using two different approaches, we then in- 
vestigate the collective excitations: we model the N phase 
both hydrodynamically and collisionlessly. The collective 
mode spectrum is discontinuous at the first-order tran- 
sition and has frequencies below the trapping frequency 
near the continuous transition. Finally, we find that the 
polarization of the N phase close to the interface has a 
strong effect on the thermal equilibration process: at low 
temperatures, energy transport through the N-SF inter- 
face is blocked. 
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II. TRAPPED SYSTEM 



p, in the N (i =|, J,) and SF phases are related by: 



A. Trapping Potential 

Experimentally, a quasi-2D regime is established by 
combining tight confinement along the z-axis weak har- 
monic trapping in the x — y-plane. The tight confine- 
ment is achieved using an one-dimensional optical lat- 
tice, which is modelled by a sinusoidal potential. The 
total trapping potential is: 



V(r) = V sin' 



where r 2 — x 2 + y 2 . If the energy cost of crossing the 
potential wells along the z-dircction is much larger than 
the Fermi energy, the inter-well tunnelling rate is negligi- 
ble and the particles are essentially confined in quasi-2D 
layers. In such a quasi-2D system a two-particle bound 
state always exists and has energy [35|, H? 



E b = 0.915^ exp (V2^£ z /aj , 



(1) 



where we have introduced to 2 — 87r 2 Vo/(mA 2 ), £ 2 — 
h/muo z and a is the 3D s-wave scattering length. Within 
mean-field theory, the bound state energy Ef, fully char- 
acterizes the interparticle interactions and can be used 
to find qualitative results for the entire BEC-BCS transi- 
tion. Therefore, as opposed to the 3D case, the crossover 
may be explored not only by tuning the scattering length 
a but also by adjusting the trapping parameters Vq and 



B. Equation of State 

Let us now introduce the mean-field formalism which 
we will use later. We follow Ref. in which homoge- 
neous systems were thoroughly studied and analytical re- 
sults were derived. This theory should be quantitatively 
correct in the BCS regime but also yields physically rel- 
evant results for the entire BEC-BCS crossover. 

In the following, we consider two homogeneous phases: 
the normal (N) and the superfluid (SF) phase. The N 
phase is taken to consist of a non-interacting mixture of 
majority (f) and minority (J.) particles. The N phase 
is said to be fully polarized (FPN) when the density of 
the minority species is zero (pi = 0) and partially polar- 
ized (PPN) otherwise. The interactions in the SF phase, 
on the other hand, are characterized by the bound state 
energy Eh of Eq. (fT]). As was found both theoretically 
(Ref. Hi) and experimentally (Refs 0, i, S fl S S 0), 
the two spin densities are equal in the SF. We use the lo- 
cal density approximation in which, at each point in the 
trap, the results for the homogeneous phases are used. 

At each point r in the trap, the pressure, P and density, 



= ^rpfCr) and PsF(r) = Z^k«. (2) 
m 2m 

Also, the densities and chemical potentials p obey the 
following relations: 

2 7 rfi 2 ( o 4 (r) irh 2 p SF (r) E b 

Hi(r) = and p S F{r) = — . (3) 

m ml 

Within the LDA in the x— y-plane, p(r) = p—muj 2 r 2 /2 
for all phases. We define the length scales Ri and 
Rsf: 



and psF(r) = 



R 



SF 



2tt#, 



(4) 



with i =1,1 and l u = ^h/mujQ. Finally, since the SF 
state is balanced, it consists of an equal | and j density 
at each point of the trap. 

The N and the SF phases are separated by an interface 
at radial position £ where two boundary conditions apply. 
First, mechanical equilibrium at the interface demands 
that the difference of the pressures at the interface must 
be compensated for by a surface-tension term [23| ■ This 
is expressed by Laplace's formula [50j : 



SF 



p i)c = ( ? 



(5) 



C 



where a is the surface tension and (•)<; means that the 
quantity is evaluated at the interface. The second bound- 
ary condition pertains to the fact that pairs of opposite- 
spin particles may cross the interface [2l|. Therefore, at 
equilibrium, the energy to add such a pair of particles in 
the N phase must be equal to the energy to break a pair 
in the SF phase, or: 



2(ms_f) c = (Ml +IH)c,- 



(6) 



In case the interface is between the SF and the FPN 
phase, no J. particles are present in the N phase. There- 
fore, in Eq. ([5]) the pressure P; can be put to zero. We 
will later see that in that case, the interface becomes im- 
permeable, that is, no particles can cross the interface, 
which is as it should be. 



C. Trap Characteristics 

Using the local properties mentioned above we can ex- 
tract global properties of the trap. We take the trap 
to contain N = iV-p + N± particles in total. The particle 
numbers of the different spin species JVj and N± are given 
by 



Q = (N 1 -N l )/(N ] +N l ) 



(7) 



3 



where Q is the population imbalance. The particle den- 
sities of Eq. (0| can be integrated to yield the particle 
numbers in the N shell, Npfi, and in the SF core, N$f'- 



a) Q=0.17 



b) Q=0.3 



N Nl = 



N SF = 



( 2 ( R SF 



em 



2#, 



(8a) 



(8b) 



with i =|, [. Also, since the SF phase is unpolarized, the 
total particle number in each spin state is given by 



N t = N Ni + Nsf/2. 



(9) 



We have now established all the ingredients necessary 
to construct all possible equilibrium states. This is done 
by solving the equations above for the variables CjRsf, 
i?l and R[. Three trap configurations turn out to be 
of relevance. First, a trap consisting of a SF core phase 
surrounded by a partially polarized normal (PPN) shell, 
which we denote by PPN+SF. A second configuration, 
denoted FPN+SF, consists of a SF core phase surrounded 
by a fully polarized normal shell [13] • Finally, the trap 
fully normal configuration is simply denoted by FN. 

In Fig. Q] the spin densities of these three trap configu- 
rations are shown. The minority spin density, pi vanishes 
in the FPN phase and pj = pi in the SF phase. 



D. Trap Characteristics without Surface Tension 




FIG. 1: The spin densities pj and p^ of the three relevant 
trap configurations as a function of the radial coordinate r. 
The density of the majority species pj is drawn using a full 
line, while the minority species py is shown with a dotted line. 
We took E^/Ef = 0.04 and varied the population imbalance 
Q. £ is the position of the interface, while _R| and Ry are 
defined in Eq. @. a) At low population imbalance (here 
Q — 0.17) we find a configuration consisting of a SF core 
surrounded by a fully polarized normal shell (FPN+SF) . b) At 
intermediate population imbalance (here Q = 0.3) we find the 
trap to consist of a SF core surrounded by a partially polarized 
normal shell (PPN+SF). c) At sufficiently high polarizations, 
the trap is fully normal (FN). 



We continue by giving some analytical results for zero 
surface tension, a = 0. The case of a nonzero surface 
tension is treated in Sect. IIII CI 

We begin with the PPN+SF configuration which, at 
the interface, is characterized by a partially polarized 
normal phase: in terms of energy this means that h < 
(Msf)j where h = (/ij — /i|)/2 is a constant through- 
out the trap. The coexistence condition Eq. ([5]) and the 
particle number equation (|8ap reduce to [36] : 



h 2 

Ajv| — Nni 



HSFEb H 
U^ ) 2 



El 



(10a) 
(10b) 



Combining these we find that R 2 SF = 2£lVN and the 
interface position £ can be obtained as a solution of the 
following equation: 

2Q 2 E F =E b (2E F -E b - m Wo 2 (f 

x (2E F -E b /2-muj 2 C 2 ) 1 (11) 

where we have defined the Fermi energy (5l| : 

E F = thj N 1/2 . (12) 

Note that all parameters of Eq. ([IT]) are global trap char- 
acteristics. 



In contrast to the PPN+SF configuration, exact an- 
alytical expressions may be obtained for the FPN+SF 
configuration. Full population imbalance at the N-SF in- 
terface implies that, at the interface, h > (^sf)^- There- 
fore, Eqs. © and (JHaJ) reduce to [3f|: 



2N 



N1 



( hsf + h 



These equations result in the following radii: 



R 2 SF = 2£lVN, 



(13a) 
(13b) 



(14a) 
(14b) 

(14c) 



and clearly £ < Rgp < R^ as it must. Note that R$f is 
independent of Q and, remarkably, all lengths are inde- 
pendent of E b . In the limit of a balanced trap (Q = 0), 

all radii equal \J 2\/Nl w ; when, on the other hand, the 
trap is fully imbalanced the interface position vanishes 
while i?| becomes equal to the Thomas-Fermi length 
Rtf=I u {&N) 1 /*. 



i? 2 = 2ll VN [1 + v/Q(V2 - 1; 
C 2 - 2£l^N ( 1 - ^/Q 
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III. PHASE DIAGRAMS 

Having established the possible trap configurations we 
present now phase diagrams of trapped imbalanced 2D 
fermion gases at zero temperature. After introducing the 
concept of total trap energy in Sect. IIII A"I we discuss the 
phase diagram without (Sect. lTHB]) and with fSect. HlTCl) 
surface tension. 



A. Total Trap Energy 

At zero temperature, whether or not a certain trap con- 
figuration is the ground state is determined by its total 
energy. Therefore, one must compare the energies of all 
three trap configurations at fixed total particle number 
N and population imbalance Q. 

By straightforward integration, we find the energy con- 
tributions of the normal shells E^i (i =t, !), the SF core 
phase E$f an d the surface tension E a to be: 



Em 

Esf 
E n 



2itC,<t. 



[C 4 + 3R 2 SF {R 2 SF 



E b N SF /2, 



Summing all the contributions, we obtain the total en- 
ergy of each of the three possible trap states. For in- 
stance, the energy of the FN configuration can be ob- 
tained by setting ( = 0, which yields: 



E,, 



Ep 
6 



((l-Q) 3 / 2 + (l + Q) 3 / 2 ). (15) 



For the PPN+SF configuration with a zero surface ten- 
sion, we find: 




Here, the interface position £ is found by solving Eq. (TlTj) . 

Lastly, £jvx must be set to zero when calculating the 
total trap energy of the FPN+SF state. We find, in the 
case of zero surface tension, 



Ep 
3 



1 



)3/2 



2-1) 



yW-1). (17) 



These expressions for the energy allow us now to obtain 
the phase diagram for zero surface tension. 



B. Phase Diagram with Zero Surface Tension 

We consider again trapped systems without surface 
tension (a = Q). In this case we are able to describe the 
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FIG. 2: Phase diagram for trapped 2D imbalanced fermion 
gases as a function of the population imbalance Q (see Eq. 0) 
and Eb/EF with a zero surface tension. At sufficiently low 
population imbalance the trap configuration (FPN+SF) con- 
sists of a SF core and a fully polarized normal shell. At in- 
termediate population imbalance, the PPN+SF trap consists 
of a SF core and a partially polarized normal shell. Finally, 
in the FN configuration, the system consists of normal par- 
ticles only. The phase boundaries Q c i, Qci and Q C 3, the 
spinodal Q C 4, and the points T and TP are analytically given 
by Eqs. (|18[l - (|23[l . Notice the logarithmic scales. 



entire phase diagram analytically. The phase diagram is 
mapped out in Fig. [2 as a function of the population im- 
balance Q and the dimensionless parameter Eb/Ep. The 
diagram is characterized by the triple point TP which, 
as we show later, has the following coordinates: 



{Eb/Ep, Q) 



V2 



TP 



(1 + V3), 1/2 



(18) 



The triple point lies at the intersection of the lines Q c \ , 
Q c2 and Q c3 , which separate the FPN+SF, PPN+SF and 
FN trap configurations. It is evident from Fig. [J] that, 
for population imbalances above 0.5, the PPN+SF is not 
stable for any interaction regime. 

The line Q C 3 cuts the line of full population imbalance 
Q = 1 at the point T which, as shown later, has coordi- 
nates: 



(E b /E F ,Q) T = (l-^, 1 



(19) 



This means that if the interaction parameter is larger 
than 1 — V2/2 then the trap always contains SF particles. 

For small enough population imbalance Q, the 
FPN+SF configuration is the one of lowest energy, 
whereas for large enough Q, the trap configuration 
is always FN. At intermediate population imbalances 
and when the interaction strength is sufficiently weak 
(Eb/Ep < 1 — \/2/2), a partially polarized shell may be 
found in the N phase. The transitions occurring along the 
lines Qd and Q c s of Fig. [JJ are discontinuous while along 
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FIG. 3: Comparison of the total trap energy Etot of three rel- 
evant trap configurations against the population imbalance Q 
in the case Ei / Ef =0.04 and with no surface tension. For low 
enough Q < Qci = 0.171 the FPN+SF configuration has the 
lowest energy. For Q cl < Q < Q c2 = 0.39 it is the PPN+SF 
configuration that is stable, while a FN trap is encountered 
in the case Q > Q C 2- The transition at Q c \ is discontinuous; 
the one at Q C 2, on the other hand, is continuous. 



Q C 2 they are continuous. In Fig. wc illustrate the na- 
ture of these transitions by plotting the total trap energy 
of the different trap configurations against the popula- 
tion imbalance Q in case E b /Ep = 0.04. For population 
imbalances below Q c \ = 0.171, the configuration of min- 
imal energy is FPN+SF while for Q c \ < Q < Q C 2 it is 
the configuration with a partially polarized normal shell 
(PPN+SF). Both configurations have metastable exten- 
sions below and above Q c i- The transition at Q C 2, on the 
other hand, is continuous and the PPN+SF configuration 
has no metastable extension above Q C 2- Approaching Q C 2 
from below, the equilibrium trap configuration is charac- 
terized by a vanishing population of SF particles. 

The dashed line in the phase diagram of Fig. [5] is 
the spinodal line of the PPN+SF configuration. It de- 
notes the line below which no metastable PPN+SF con- 
figurations exist. Qc4 is an exceptional spinodal in the 
sense that, approaching the Q C 4 line, the metastable 
PPN+SF configurations have a vanishing population of 
minority atoms in the N phase; therefore, on the Qd 
line, the PPN+SF configurations have essentially become 
FPN+SF. Note also that there are no spinodal lines for 
the FN or FPN+SF configurations. 

For all lines in Fig. [5] we have obtained analytical ex- 
pressions. Firstly, the critical line Q C 2 can be easily found 
since we know that £ vanishes at the transition. Putting 
C = in Eq. (HU) yields [H : 

Secondly, the phase boundary between FPN+SF and N 



FIG. 4: The different radii (, and Rsf in units of 

Rtf = ^(S-ZV) 1 / 4 as a function of the population imbalance 
Q in case Eb/EF =0.04 and with a zero surface tension. 

is readily found by equating Eqs. (fT5|) and (fTT)) : 

Eb j(l + Q% 2 (V2-l))-(l-Q c3 )^-(l + Q c3 )^ 
~E~p ~ 3(1 -Q c3 ) 

(21) 

Now, the triple point TP and the point T follow from 
Eqs. (|20p and (f2Tj) . The spinodal line Q C 4 denotes the line 
where, at the interface, the FNP becomes partially polar- 
ized and thus h = (^sf)^ and with Eqs. (|10ap and ()13a|) . 
we arrive at [381 ]: 

<H^)(t) 2 - <22) 

Finally, wc find (numerically) that Q c \ cx {E b /Ep) 2 . 
Given the fact that Q c \ goes through the triple point 
TP, we conjecture that the phase boundary is given by: 

Qcl = [2-^(l + V3)P (SO • (23) 

Eq. (|23p agrees very well with the numerically obtained 
results. 

The phase diagram of Fig. [2] differs from that of Fig. 4 
in Ref. HH. Indeed, the line shown there as the phase 
boundary between the FPN+SF and the PPN+SF in 
that paper corresponds to line Qd in our Fig. [21 A 
comparison of the energies of the possible configurations, 
however, shows Q C 4 to be a spinodal line and Q c \ to be 
the phase boundary instead. In such cases, when more 
than one trapped density configurations satisfy the num- 
ber and imbalance constraints, comparing their energies 
to determine the true ground state is crucial. 

The changes to the trap configurations at the two tran- 
sitions Qd and Q C 2 are illustrated in Fig. [4] where we plot 
the different radii as a function of the population imbal- 
ance Q, again for E b /Ep = 0.04. For low enough popu- 
lation imbalance, Q < Q c i, the interface of the FPN+SF 
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configuration moves inwards with increasing population 
imbalance; the radii vary according to Eq. fT4")) . Note 
that the radius i?| is not defined for the FPN+SF config- 
urations as the N phase is fully polarized. The first-order 
transition at Q c \ introduces discontinuities in the radii. 
On further increasing the population imbalance above 
Qd, the radius at which the N-SF interface lies, £, de- 
creases and finally vanishes at the continuous transition 
Q C 2, during which all radii vary continuously. 

C. Phase Diagram with Surface Tension 

We have also numerically calculated the phase dia- 
grams for 2D trapped imbalanced fermion with nonzero 
N-SF interface tension (a ^ in Eq. ([5])). No exact ex- 
pression for the surface tension is available at this point. 
Nevertheless, dimensional arguments at unitarity and ex- 
plicit calculation in the BEC limit [52j yield 

o-(r)=r,H 2 p 3 s / *(r)/m. (24) 

Here r\ is an arbitrary dimensionless constant, which we 
take equal to unity. While the analysis leading to Eq. (|24|) 
is not strictly valid away from unitarity, we believe that 
the qualitative features of the phase diagrams are not 
altered by these corrections. 

We present our phase diagram for a fixed particle num- 
ber N = 10 4 in Fig. [5] As is expected and clearly seen by 
comparing Figs. [2] and [3 the introduction of surface ten- 
sion strongly suppresses the appearance of the SF phase. 
Specifically, the triple point TP of Fig. [5] has moved to 
vanishingly small values of the two variables Ef, fEp and 
Q, removing the PPN+SF configuration from the dia- 
gram. We have also found the TP to be at vanishingly 
small values for the case N = 10 8 . 

The effect of a nonzero surface tension is therefore to 
enlarge the region where the FN configuration is present. 
As was also found for the case of zero surface tension, 
the phase boundary Q c % is discontinuous. Note that the 
phase boundary Q c ^ does not seem to intersect the line 
of full polarization Q — 1, so that the point T of Fig. [2] 
is not present at finite surface tension. 

IV. COLLECTIVE EXCITATIONS 

In this section we study the collective excitations of the 
imbalanced 2D fermion gases. The analogous excitations 
of a 3D system were already studied in the experiments 
of Ref. |24j. We stress the role of the continuous and dis- 
continuous transitions as presented in last section. First, 
in Sect. IIV Al we take the approximation that both the 
N and the SF phases in the trap behave hydrodynami- 
cally. This allows us to treat a permeable interface, that 
is, an interface at which particles may go over from the 
N to the SF phase, as is present in case of a PPN+SF 
configuration. Then, in Sect. IIVB| we treat a collision- 
less N phase and are therefore able to treat exactly the 




E b/ E F 

FIG. 5: Phase diagram for trapped, 2D imbalanced fermion 
gases, as a function of the population imbalance Q (see 
Eq. 0) and E^/Ef with a nonzero surface tension and 
N = 10 4 . At sufficiently low population imbalance the trap 
configuration (FPN+SF) consists of a SF core and a fully po- 
larized normal shell while the FN configuration (dark grey 
array) consist of all normal particles. The Q c z phase bound- 
ary is of first order. Notice the logarithmic scales. 

collective excitations of the FN and, also of the FPN+SF 
configurations. 

A. Hydrodynamic Phases 

We begin by calculating the expressions for the bulk 
density fluctuations in the trap for both the SF and the 
N phase, assuming their dynamics can be described by 
hydrodynamics. Since we are interested in the lowest- 
energy frequencies of the trap when subjected to a weak 
excitation, we consider only the linearized hydrodynamic 
equations as usual. 

Experimentally realized ultracold Fermi gases are 
highly compressible and, at each position, their densi- 
ties deviate from their equilibrium value by an amount 
Spj with j =T)J. and SF. Henceforth, we will denote 
the deviations from equilibrium of the chemical poten- 
tials and the pressures by dp and 5P respectively. An 
external perturbation of the system will induce the sys- 
tem to move at (position-dependent) velocity v in such a 
way as to sat isfy the linearized Euler and the continuity 
equations H E& 11 HI E| : 

mpjdtVj = -VSP 3 - 6pW, (25a) 
d t 5 Pj = -V • ( ft v), (25b) 

where we only allow the particles to move in the x-y plane 
such that we can consider V(r) — muj 2 r 2 /2. Combining 
Eqs. (|25a|) and (|25b|) and writing Spj(t) -> Spje~ iuJt , it 
follows that: 

mu! 2 Spj = - mV 2 [c 2 8p.j\ - V • (SpjW), 

where c is the position-dependent velocity of sound which 
is introduced through the use of the Gibbs-Duhem rela- 
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tion [§, [l(| ■ We continue by performing the transforma- 
tion Spj(r,9) — > 5pj(r)e lte where £ is a positive integer 
or zero, so as to get: 



= [l-rf] 



Sf 2 



^1 



2r 



dr~ 



+ 2uj 2 5pj/u>Q 



with r = r/Rj. Solutions to this equation can be 
straightforwardly obtained in terms of the hypergeomet- 
ric function F; moreover, they must be regular at r = 
when j = SF and at r = Rj when j =|, J,. Moreover, 
from Eq. ([3]) , it is clear that density deviations should be 
proportional to the chemical potential perturbations, or 
Sp oc Spj. Thus 



SpsF oc r F 
Spi oc r £ F 



a + ,a-,a»,(r/Rs F y) 



a 



(26a) 
(26b) 



with i =t, |, a = £+1 and 2a± = e+l±[£ 2 +2uj 2 / 'w 2 ] 1 / 2 . 

As in Sect. Ill Bl we must supplement these bulk equa- 
tions by boundary equations at the N-SF interface. Me- 
chanical equilibrium, as was expressed in Eq. for the 
equilibrium situation, is also demanded to be valid away 
from it. This results in the following expression: 



PsfSpsf - Pt^T _ Pi 5 Pi + s Cd r (Psf - P\ - P[) 



= <j5( 



£ 2 -l 



+ 



SmadpsF 



(27) 



C 2 2-K PSF £i) 2kV Psf C 

As mentioned before, particles may cross the permeable 
interface. At equilibrium, this was encoded by the chemi- 
cal potential balance Eq. ((5]); away from equilibrium, this 
gives: 

<^T + S H ~ ZSpsF = 5( d r {2p SF - /x T - p,±). (28) 

Lastly, the interconversion of particles between the N and 
the SF phase is not only restricted by chemical equilib- 
rium (Eq. (f2"8")) ). but also the continuity equation must 
be satisfied at the interface. Therefore, the fluctuating 
interface position S( depends on the velocities in the N 
and the SF as follows [421 ]: 

e c • (2vipi - v S fPsf) ( = (2pi - p S F)(d t 5(, (29) 

with i =|, I and is a unit vector perpendicular to the 
interface and directed towards the N side. Assume now 
that the N phase is fully polarized, such that (pi)^ = 0. 
It is then easily shown that ~vsf = Vf = dtSC, implying 
that the interface is impermeable; the interface cannot 
move by interconversion of particles between the N and 
SF phase. In what follows we show that the collective 
excitations are strongly affected by the permeability of 
the interface. Taking the time derivative of Eq. ([29]) and 
using the Euler equation yields: 

mu?(2pi - psf) S( = 2p { d r {Spi) - pspd r {Spsf) ■ 

Prior to presenting our results, we outline the theory 
for treating a collisionless N phase. 



B. Collisionless Normal Phase 

We briefly discuss how to treat the collisionless N 
phase in a harmonic trap. We shall use the formalism 
presented in Refs. [23| and|43| and refer to these articles for 
details. In what follows, we treat the N phase as a non- 
interacting gas; this is valid when the N phase is fully 
polarized at the interface, that is, for all the FPN+SF 
configurations. As mentioned before, Eq. (|2"5|) then im- 
plies that the interface is impermeable to particles. 

The SF core, on the other hand, still behaves hydro- 
dynamically, so that the chemical potential perturbation 
still satisfies Eq. (|26a|) . On the other hand, Pauli block- 
ing in the N gas strongly suppresses interactions, leading 
to long local equilibration times and rendering a hydro- 
dynamic description inapplicable. We must instead solve 
the Botlzmann-Vlasov equation for the full distribution 
function [23| . The hydrodynamic description uses four 
numbers at each point, p(v,t) and v(r, t) while the de- 
scription of the N gas is based on the full distribution 
function /(r,v,i); matching the two at the interface is 
nontrivial. The problem has been solved for the three- 
dimensional system in Ref. [23j . 

In brief, it is necessary to impose two boundary con- 
ditions: First, impermeability means that the interface 
velocity must equal the SF velocity: 



d r (5psF) =muj 2 5(. (30) 



Secondly, in the condition for mechanical equilibrium, 
Eq. (j27|) . we must set P± = pi = Q and replace the pres- 
sure term P|- by the momentum flux tensor <5nJ r . In 2D 
and for radially symmetric perturbations, a generaliza- 
tion of the calculations of Refs. 23, 43 leads to: 



tkj5((R 2 



2\3/2 



C 2 ) 



rrj" 



el* 2 



1 rf x 3 cot(wr/2) 



with 



arctan 



2 X ( ^(P T /C) 2 - i - iA\Mt/C) 2 



Ab- 



solving these equations with respect to u> gives the fre- 
quencies of the system's breathing modes. 



C. Results 

In Figs. [5] and [7] we present the spectra for the collec- 
tive excitations of trapped imbalanced fermion gases in 
2D taken into account a hydrodynamic (Sect. IIV A[) and 
collisionless N phase (Sect. IIV Bp . respectively. We be- 
gin by discussing the fully hydrodynamic case of Fig. [5] 
which shows the results for a system with a zero sur- 
face tension (er = 0) and Eb/Ep = 0.04. We dis- 
tinguish two types of modes: the in-phase (IP) modes 
(full lines) and the out-of-phase (OOP) modes (dashed 
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FIG. 6: The collective mode frequencies as a function of the 
population imbalance Q in case the trap is fully hydrody- 
namic, E^/Ef = 0.04 and no surface tension is present. We 
show the in-phase (full lines) and out-of-phase (dashed lines) 
spectra of the £ = 0, 1 and 2 modes. For low population im- 
balance Q < Qd the system has a FPN+SF configuration 
while for Q c \ < Q < Q C 2, it is the PPN+SF configuration 
and it is FN for Q > Q C 2- A fully hydrodynamic formalism 
was used in order to be able to incorporate the permeable 
interface for the PPN+SF configuration. 



line) [23J, |44|. The OOP modes are unique to a two- 
component system. In-phase and out-of-phase refer to 
the relative motion of the N-SF interface with respect to 
the outer boundary of the N shell (see Fig. 1 in Ref. l23l). 
Also, we can also distinguish between three regimes de- 
pending on the trap configuration: from the phase dia- 
gram of Fig. [21 it follows that for Q < Q c i = 0.17, the 
trap configuration is FPN+SF, while it is PPN+SF for 
Qd < Q < Qc2 = 0.39 and FN for Q > Q c2 . At zero 
and full population imbalance, the single-species hydro- 



dynamic spectrum cu 



2n(n 



1)] is recov- 



ered. For small but finite values of Q, OOP modes with 
low frequencies appear (note that we only show t — 1 
and 1 = 2). As is evident in Fig. discontinuities 
in the trap-characteristic radii appear at the first-order 
transition Q = Q c i\ the presence of discontinuities in 
the collective mode spectrum is therefore not surprising. 
More surprising is the appearance of additional spectral 
lines for values Q c \ < Q < Q C 2- These extra modes are 
due to the presence of spin J, particles in the N shell, 
which gives rise to an extra degree of freedom. Since the 
transition at Q c i is continuous, no discontinuities appear 
in the spectrum there at that point. For Q > Q C 2 the 
trap is fully N and the spectrum becomes independent of 
the imbalance. Finally, notice that the breathing mode 
(black line) which has a frequency of exactly u> = 2luq 
when Q < Qci, becomes lower at Q c \ and vanishes at 
Q C 2i that is when the SF core phase disappears. 

In Fig. [7] we consider the lowest two breathing modes of 
a system with nonzero surface tension, with N = 10 and 
with interaction parameter E\,jEp = 0.5; we treat the 
N particles collisionlessly. At zero polarization, we find 



CO 4 
CO 






I 

Damped 



(2,0) 



(1,0) 



0.5 

Polarization Q 



■C3 



1.0 



FIG. 7: The two lowest breathing mode frequencies as a func- 
tion of the population imbalance Q in case Et/EF = 0.5 and 
with a collisionless normal trap with a nonzero surface ten- 
sion. For population imbalance below Q C 3 = 0.81 the system 
has a FPN+SF configuration while it is FN when Q > Q C 3. 
The dashed region indicates collisionless interfacial damping. 



the spectrum of the two lowest hydrodynamic breath- 
ing modes 2u>o and 2\/3wo, while at full trap imbalance, 
we recover the collisionless breathing mode frequencies 
2ujo and Au>o. As was also found in Fig. [HI the lowest 
breathing mode of the FPN+SF configuration in Fig. [7| 
lies exactly at 2tuo for all Q; we consider this feature to 
be an indication of the correctness of our formalism, even 
in the presence of a nonzero surface tension. For popula- 
tion imbalances Q such that 0.17 < Q < 0.27, the second 
mode is damped. This damping is caused by fermions 
which are resonantly driven at the interface as explained 
at length in Refs.[23| andl43l Note also the weak discon- 
tinuity at the Qc3 occurs when the trapped configuration 
changes from FPN+SF to FN. The amplitude of this dis- 
continuity changes with the interaction parameter E^. 



V. ENERGY TRANSPORT 

We discuss now how heat is transferred throughout a 
trapped population- imbalanced fermion gas. Specifically, 
we focus on the heat transport through the N-SF inter- 
face [Til [2l|. We find that at low temperatures, energy 
flux is exponentially suppressed due to the presence of a 
SF gap. This due to the suppression of particles crossing 
the interface at low temperatures, since heat is transmit- 
ted through the N-SF interface when thermally excited 
particles penetrate the interface. 

To understand this, consider a N particle incident on 
to the interface. 

If the N phase is partially polarized (for PPN+SF 
configurations), there are four scattering possibilities: 
Andreev reflection, specular reflection and hole-like and 
particle- like transmission [46| . A particle undergoing An- 
dreev reflection pairs up with a particle of opposite spin 
to form a Cooper pair in the SF and therefore leaves be- 
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hind a hole in the N phase. Particles penetrating the SF 
must have a minimum energy A^ — h, where is the 
SF gap at the interface [HJ. However, at sufficiently low 
temperature T, this is rare because of the low statistical 
weight e ~( A c-' l )/ fcsT < 1. 

If, on the other hand, the N phase is fully polarized 
at the interface (for FPN+SF configurations), both An- 
dreev reflection and quasiparticle transmission are sup- 
pressed so that all N particles must be specularly re- 
flected off the interface. Therefore no energy transport 
through the interface is possible for FPN+SF configura- 
tions [531 ] . 

We now quanitfy this suppression of energy transport 
in the case of a PPN+SF configuration. To do this we 
calculate the heat resistivity of a N-SF interface on the 
basis of the Bog oliubov-de- Gennes equations, referring to 
Refs . EH! and l45l for details. A brief outline of the followed 
method is now given. We model the N-SF interface by 
its step-like behavior of the gap function: in the N phase 
A = 0, while according to mean-field theory and coexis- 
tence condition (|10ap [3a |: 

A c = y/2h* + E*/2, 

in the SF. Note that since we consider PPN+SF config- 
urations, h < Et,(l + \/2)/2. Thermally excited particles 
and holes which are incoming on the N-SF interface are 
described by quasiparticle wave functions. Matching the 
wave functions and their derivatives at the interface al- 
lows us to calculate the transmission coefficients. We 
then calculate the heat conductivity K of the N-SF in- 
terface which is the net heat flux per temperature differ- 
ence across the interface. We find that in case a small 
temperature bias between N and SF phase exists, the 
thermal equilibration process will be blocked by a ther- 
mal conductivity that decays exponentially fast as the 
temperature falls. In the limit of very weak interactions 
Eb/nsF "C 1 and at low temperatures ksT <§C A^ — h, 
one can use the Andreev approximation to obtain the 
following analytical expression for the two-dimensional 
case dalU El: 

^_ _ e-( A <- gh >/ fc * T V?(A c -(7/i) 2 
~i^~£^ ± k B Ty/2k B TAs ' 1} 

It is clear that the dominant factor at small temperatures 
is e~( A< ~ h ^ kBT . We have found that the approximation 
of Eq. pip is in very good agreement with numerically 
obtained results, even for large values of Eb/ hsf- 

We therefore conclude that, at sufficiently low temper- 
atures, thermal equilibration between N and SF phases at 
different temperatures will be considerably slowed down 
due to the SF gap. 

VI. DISCUSSION 

Although our formalism is only quantitatively correct 
in the BEC regime, we expect that it also yields qualita- 



tively realistic information throughout the entire BEC- 
BCS cross-over 49]. The condition Eb <C ^sf is gener- 
ally correct in trap center. 

We did not take into account the posibility of Fuldc- 
Ferrel-Larkin-Ovchinnikov (FFLO) states [H, In 
Ref. [33[ it was found that such FFLO phase would ap- 
pear between the PPN and the SF. Our results therefore 
lead us to speculate that the FFLO phase would be sup- 
pressed by the presence of a nonzero surface tension and 
the collective excitations as presented in Fig. would still 
be correct. 

In future work, we would like to investigate the finite- 
temperature modifications to the results obtained here. 
At nonzero T, the long-range order in the SF is absent in 
2D. This implies the presence of a Berezinskii-Kosterlitz- 
Thouless transition HO, HH [HJ . 



VII. CONCLUSION 

We have studied population imbalanced Fermi mix- 
tures under quasi two-dimensional confinement. Us- 
ing results previously published for the homogeneous 
case and using the local density approximation, we have 
shown that there are three relevant configurations of the 
trapped cloud: fully normal (FN), a partially polarized 
normal cloud surrounding a superfluid core (PPN+SF) 
and a fully polarized normal cloud surrounding a super- 
fluid core (FPN+SF). 

We first show that, in the absence of surface tension, 
the system is completely specified by two numbers: the 
polarization, Q = (iVf — N{)/(N-[ + ATjJ (where N. t the 
total number of particles of species i), and the interaction 
parameter Eb/Ep (see Eqs. ([TJ and (|12p for definitions). 
Thus, the phase diagram we obtain (Fig. [2]) is universal, 
in the sense that it only depends on a small number of 
parameters. This phase diagram has the property that, 
for E b /E F > 1 - \/2/2 (see Eq. (fT9"]0 there is no critical 
Q c above which the SF disappears; in other words, above 
this Eb/Ep, there is SF in the trap for all polarizations. 

A finite surface tension destroys this universality, so 
we have studied systems containing specific numbers of 
particles. The qualitative changes in the system's ground 
state are clear: the SF and PPN phases are suppressed. 

We then turn to the calculation of collective excitation 
modes of the system. We study these in two approxima- 
tions. 

First we take both the SF and the N parts to be- 
have hydrodynamically. In the FPN+SF case we find 
out-of-phase modes, typical of two-component systems 
(and which we have also found in previous work on 3- 
dimensional systems). In the PPN+SF case, on the other 
hand, the presence of three components (SF, PPN and 
FPN) results in an extra set of modes. 

We then study the more physically relevant case where 
the fully polarized N gas displays collisionless behaviour. 
In this approximation we calculate the collective mode 
frequencies and find the expected limits, as well as damp- 



10 



ing (related to Landau damping) for a range of polariza- 
tions. 

Finally we discuss the effect of the interface on ther- 
mal transport through the cloud. This would be im- 
portant in studying how the trapped cloud equilibrates 
starting from an initially nonuniform temperature distri- 
bution, such as might be obtained while the system is 
being cooled. Based on previously published work, we 
find that thermal transport is strongly suppressed at low 
temperatures due to the presence of the gap in the exci- 
tation spectrum inside the SF. 
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